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ABSTRACT 

Making use of scaling relations between the central and the total galaxy luminosity of 
a dark matter halo as a function of the halo mass, and the scatter in these relations, 
we present an empirical model to describe the luminosity function (LF) of galaxies. 
We extend this model to describe relative statistics of early-type, or red, and late- 
type, or blue, galaxies, with the fraction of early type galaxies at halo centers, relative 
to the total sample, determined only by the halo mass and the same fraction in the 
case of satellites is taken to be dependent on both the halo mass and the satellite 
galaxy luminosity. This simple model describes the conditional luminosity functions, 
LF of galaxies as a function of the halo mass, measured with the 2dF galaxy group 
catalog from cluster to group mass scales. Given the observational measurements of 
the LF as a function of the environment using 2dF, with environment defined by the 
galaxy overdensity measured over a given volume, we extend our model to describe 
environmental luminosity functions. Using 2dF measurements, we extract information 
related to conditional mass function for halos from extreme voids to dense regions in 
terms of the galaxy overdensity. We also calculate the probability distribution function 
of halo mass, as a function of the galaxy overdensity, and use these probabilities to 
address preferred environments of red and blue galaxies. Our model also allow us 
to make predictions, for example, galaxy bias as a function of the galaxy type and 
luminosity, the void mass function, and the average galaxy luminosity as a function 
of the density environment. The extension of the halo model to construct conditional 
and environmental luminosity function of galaxies is a powerful approach in the era of 
wide-field large scale structure surveys given the ability to extract information beyond 
the average luminosity function. 

Key words: large scale structure — cosmology: observations — cosmology: theory — 
galaxies: clusters: general — galaxies: formation — galaxies: fundamental parameters 
— intergalactic medium — galaxies: halos — methods: statistical 



1 INTRODUCTION 

An important aspect of understanding underlying astro- 
physical reasons for galaxy formation and evolution involves 
studying the relative distribution of red (early-type) and 
blue (late-type) galaxies, as a function of the galaxy environ- 
ment. Observational measurements of this so-called density- 
morphology relation suggest evidence that early-type galax- 
ies are predominantly found in dense environments such as 
galaxy groups and clusters (Dressier 1980; Goto et al. 2003) 
Evidence also suggests that the formation of early-type 
galaxies predates the formation of galaxy clusters (Dressier 
et al. 1997). In the era of wide-field galaxy surveys such as 



the 2dF Galaxy Redshift Survey (2dFGRS; CoUess et al. 
2001) or the Sloan Digital Sky Survey (SDSS; York et al. 
2000), detailed statistics on galaxy types and their environ- 
ments allow one to now construct detailed models related 
statistics of galaxy types and where these galaxies are lo- 
cated. These models may then aid in explaining the under- 
lying reasons for the occurrence of galaxy types and their 
preferred environments with initial conditions given by the 
primordial density fiuctuations and cosmological parameters 
that determine the expansion. 

Following these lines, numerical and semi-analytical 
models of galaxy formation are generally pursued to model 
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and understand galaxy statistics including luminosity func- 
tions, occurrence of galaxy types, their spatial distribution, 
and clustering properties. While initial conditions and back- 
ground cosmology are known adequately, these techniques 
are yet to describe statistical measurements of the galaxy 
distribution with reasonable accuracy. The main unknown 
here comes from our limited understanding of gastrophysics 
involving how baryons cool to form stars and how the for- 
mation of stars and their evolution, including the stellar end 
products, may lead to various feedback processes that affect 
subsequent starformation. The standard text book descrip- 
tion of galaxy formation involves gas first heating to virial 
temperature during the formation of dark matter lialos and 
subsequently cooling at halo centers to form stars (Rccs & 
Ostriker 1977; White & Rees 1978) A characteristic scale in 
galaxy formation is then related to the amount of gas, ba- 
sically within the "cooling radius" , that can cool within the 
Hubble time given a mass scale for the halo. Semi-analytic 
models of galaxy formation (e.g., Benson et al. 2001) and di- 
rect hydrodynamical simulations of the galaxy distribution 
(e.g., Kay et al. 2002) generally over predict the number of 
galaxies both at the low-end and the high-end of the galaxy 
luminosity function. The bright-end of the LF is always as- 
sociated with the "over cooling problem" in numerical sim- 
ulations (Balogh et al. 2001), where hot gas cools rapidly to 
form luminous galaxies at halo centers. The faint-end prob- 
lem, involving a lack of faint galaxies in the observed LF, is 
generally explained as due to a feed back process during the 
era of reionization (Barkana & Loeb 2001) when primordial 
galaxies started to form. Models are generally evoked to ex- 
pel gas from halos, such as through heating associated with 
reionization or a first generation of supernovac (e.g.. Bullock 
et al. 2000; Benson ct al. 2002), but the gas expelled from 
small halos settle eventually in more massive halos and cool 
to form bright central galaxies with luminosities exceeding 
those observed (e.g., Benson et al. 2003). 

Ignoring galaxy growth through continuous cooling of 
hot gas, in Cooray & Milosavljcvic (2005a), the central 
galaxy luminosity growth, as a function of the halo mass 
was explained based on a simple description for galaxy merg- 
ing to the halo center with an efficiency determined by the 
dynamical friction alone. Those models also load to a char- 
acteristic scale in the galaxy formation reflected in terms of 
a flattering of the relation between central galaxy luminos- 
ity and halo mass, or Lc{M) relation, as the halo mass is 
increased. This characteristic luminosity is associated with 
a critical mass scale when the dynamical friction time scale 
becomes close to or exceed the Hubble time. In Cooray & 
Milosavljevi'c (2005b), a model for the LF of galaxies that 
relied primarily on this relation, and to a lesser extent on the 
relation between total galaxy luminosity in a given halo and 
it's halo mass, was used to show that this characteristic lumi- 
nosity is same as L* in the LF, when the LF is described with 
the Schechter (1976) form of $(L) oc (L/L*)" exp(-L/L*). 
The success in describing the LF of galaxies using the Lc{M) 
relation and it's scatter led to the conclusion that the L* in 
the luminosity function is not a reflection of efliciency as- 
sociated with gas cooling as has been argued in the past 
based on traditional models of galaxy formation dominated 



by hot gas cooling in dark matter halos (e.g., Dekel 2004). 
The characteristic luminosity is rather due to decreasing ef- 
ficiency of dissipationless merging of galaxies to a central 
galaxy as hierarchical structure formation builds up mas- 
sive parent halos. 

Other evidences for a departure from traditional ideas of 
galaxy formation and evolution comes from numerical sim- 
ulations, where some simulations now suggest that gas, as 
dark matter halos virialize, never heat to the "virial" tem- 
perature completely, but rather, shock heating during virial- 
ization forms a bimodal temperature distribution (Dckol & 
Birnboim 2004; Keres et al. 2004). Binney (2004) suggested 
that only the colder component cools to form a galaxy, while 
the hotter component remains at the same temperature. 
There is a lower characteristic scale in galaxy formation as- 
sociated with the mass scale where most gas is never shock 
heats and remains at the virial temperature, but rapidly 
cools at the center to form a galaxy. One-dimensional numer- 
ical simulation suggests this lower mass scale is (1 — 6) x 10^^ 
Mq (Dekel 2004). In this description for galaxy formation, 
galaxies in more massive halos can only grow in luminosity 
only through mergers with other galaxies. The dynamical 
friction process involved with merging produces a consis- 
tent Lc(M) relation that agrees with observations (Cooray 
& Milosavljcvic 2005a). The same relation, when combined 
with the mass function, leads to the LF, and can be fitted 
with the Schechter (1976) form The exponential drop-off of 
the LF at the bright-end is reflection of the scatter in the 
Lc(M) relation (Cooray & Milosavljevic 2005b). 

Hero, we extend the model of Cooray & Milosavljevic 
(2005b) to describe galaxy statistics measured by the 2dF- 
GRS survey. The main advantage of using 2dFGRS data is 
the availability of LF measured as a function of the galaxy 
type, and the environment, measured in terms of the galaxy 
ovordcnsity over the volume determined by size scale of 8 
/r^ Mpc (Croton et al. 2004). The 2dFGRS data also aUow 
measurements of the conditional luminosity function (CLF; 
Yang et al. 2003b), the luminosity function of galaxies as a 
function of the halo mass (Yang ct al, 2005). Our models on 
the CLF can be directly compared to these measurements 
and interesting information on the relative distribution of 
galaxy types can be extracted from the data. 

While not complicated as semi-analytical models of 
galaxy formation, the empirical modeling approach utilized 
here has the advantage that one is able to understand main 
ingredients that shape the CLFs easily. The approach builds 
upon attempts by Yang et al. (2003b) to describe the LF 
with CLFs, but assuming Schechter forms a priori for the 
CLF, and approaches that are built upon the halo model 
for the galaxy distribution (Cooray & Sheth 2002), but now 
extended to discuss conditional functions (Zheng et al. 2004 
with the stellar mass function and Zehavi et al. 2004 in the 
case of CLFs). While the halo model has been successfully 
used to describe statistics of the dark matter field (such 
as clustering: Seljak 2000, Peacock & Smith 2000 or weak 
lensing: Cooray et al. 2000, Cooray & Hu 2001), and ba- 
sic properties of galaxy clustering (e.g., Scoccimarro et al. 
2001), it is useful to consider more applications of this tech- 
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nique which can provide information on underlying physics 
related to the galaxy distribution. 

Here, we will separate our discussion to central and 
satellite galaxies and make few assumptions as possible 
through simple model descriptions. The motivation for the 
separation of galaxies to these two divisions are numerous: 
from the theoretical side, a better description of the galaxy 
occupation statistics is obtained when one separates to cen- 
tral and satellite galaxies (Kravtsov et al. 2004) , while from 
observations, central and satellites galaxies are known to 
show different properties, such as color and luminosity (e.g., 
Berlind et al. 2004). Our goal here is to consider an ana- 
lytic description for the LF with built in model ingredients 
that recover the observations. We then argue that instead of 
attempting to understand mass-averaged statistics such as 
the LF, it may be best to reproduce main ingredients that 
shape CLFs with numerical and semi-analytical models of 
galaxy formation in order to understand underlying physics. 
In the case of CLFs, we will argue that the main ingredient 
is the Lc(M) relation, and in the case of galaxy types, a 
model on the fraction of early-type and late-type galaxies 
as a function of the halo mass. If these relations and their 
observed scatter can be explained with simple physics, then 
it is guaranteed that the LF would be recovered. 

To model the environmental LF (Croton et al. 2004), we 
need information on the mass function of dark matter halos 
that corresponds to the environment of interest, whether it is 
a void or a dense region. Since this information is not directly 
available under a simple model (see, for example, Mo et al. 
2004 and the approach there that utilized numerical sim- 
ulations), here we make use of the observed measurements 
to extract information on these conditional mass functions. 
These mass functions, as well as related probabilities, then 
provide us with general information on how blue and red 
galaxies are distributed in the Universe. 

The paper is organized as follows: In the next section, 
we will outline the basic ingredients in the empirical model 
for CLFs and present a comparison to measurements using 
2dF data (from Yang et al. 2005), both in terms of aver- 
age number of galaxies and in terms of galaxy types. In 
Section 3, we will describe the LF and in Section 4 the envi- 
ronmental LF based on the conditional mass function. Using 
data from Croton et al. (2004), we will extract information 
on the conditional mass function and various statistical mea- 
surements related to the relative distribution of early- and 
late- type galaxies. We will also study the void mass function 
and compare with predictions in the literature. We conclude 
with a brief discussion of our main results in § 5. Through- 
out the paper, we assume the concordance cosmological pa- 
rameters consistent with WMAP data (Spergel et al. 2003). 
Throughout this paper, to be consistent with observations, 
we take the scaled Hubble constant to be /i = 1, in units of 
100 km s"^ Mpc"^ 



2 CONDITIONAL LUMINOSITY FUNCTION: 
EMPIRICAL MODEL 

In order to construct the luminosity function (LF), we follow 
Cooray & Milosavljevic (2005b). The conditional luminosity 
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Figure 1. Various galaxy fractions (central, satellites, late-types 
and early-types) as a function of the halo mass. The central galaxy 
fraction falls below 1 at the low mass end as we account for the 
fact that galaxy formation is inefficient at this low end with the 
function fc{M) such that not all dark matter halos host a galaxy 
at the low-end of the mass function. The early- and late-type frac- 
tion of satellite galaxies also dependent on the galaxy luminosity. 
Here, we show the fractions for satellites with luminosity of 10^'' 



function (CLF), denoted by ${L\M), is the average number 
of galaxies with luminosities between L and L + dL that 
reside in halos of mass M (Yang et al. 2003b). First, we 
separate the CLF into terms associated with central and 
satellite galaxies, such that 



<I>(L|M) = 

ML\M) = 



^c{L\M)+^s(L\M) 

UM) logio [L/Le(M)]2 



exp 



V2^1n(10)EL 
^(M)L^(*^V.(i) 



2S2 



(1) 



Here Lc{M) is the relation between central galaxy lumi- 
nosity of a given dark matter halo and it's halo mass, while 
ln(10)E is the dispersion in this relation. The central galaxy 
CLF takes a log-normal form, while the satellite galaxy CLF 
takes a power-law form in luminosity. Such a separation de- 
scribes the LF best, with an overall better fit to the data in 
the K-band as explored by Cooray & Milosavljevic (2005b). 
While a previous attempt to describe CLFs, as appropri- 
ate for 2dFGRS, involved a priori assumed Schechter (1976) 
forms, we believe the description here is more appropriate. 
Our motivation for log-normal distributions also come from 
measured conditional LFs, such as galaxy cluster LFs in- 
cluding bright galaxies, where data do require an additional 
log- normal component in addition to the Schechter (1976) 
form (Trentham & TuUy 2002). Similarly, the stellar mass 
function, as a function of halos mass in semi-analytical mod- 
els, is best described with a log-normal component for the 
central galaxies (Zheng et al. 2004). 

In the next few subsections, we will describe in detail 
other parameters associated with the CLF and how numer- 
ical values for these parameters are obtained. First, we dis- 
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Figure 2. Central galaxy luminosity as a function of the halo 
mass as appropriate for 2dFGRS fej-band. The data points at the 
low mass end are from Vale & Ostriker (2004; filled diamonds) 
while at the high end are from Yang et al. (2005). The solid curve 
is the relation obtained in Vale & Ostriker (2004) by converting 
the 2dFGRS luminosity function to extract the plotted relation 
based on the sub-halo mass function. We use this relation to con- 
struct our the CLF, and include a scatter in this relation in our 
description of central galaxy CLFs. The dashed-line is the total 
luminosity in a given halo. This relation is needed to determine 
the CLF of satellites. For comparison, we also show the central 
and total-luminosity as a function of halo mass in the K-band, as 
obtained by Cooray & Milosavljevic (2005b) based on a variety 
of observational measurements on the central and total galaxy 
luminosity as a function of halo mass (e.g., with weak lensing: 
Yang et al. 2003a; direct measurements: Lin & Mohr 2004, Lin et 
al. 2004). The Yang et al. (2005) data are direct measurements 
of central luminosities based on the 2dFGRS galaxy group cata- 
log. The remarkable agreement between these measurements and 
the relation extracted by Vale & Ostriker (2004) suggests that the 
Lc{M) relation, as appropriate for the 6j-band is well established 
at halo mass scales above ~ lO^'^ -^0- 



cuss the CLF of central galaxies and then move on to dis- 
cuss satellites. We will end this section with a comparison 
to CLFs measured in 2dFGRS by Yang et al. (2005). 



2.1 Central Galaxies 

In our description for CLFs (Eq. 1), central galaxies have 
a log-normal distribution in luminosity with a mean deter- 
mined by the Lc{M) relation. The scatter in the Lc{M) rela- 
tion is captured through the dispersion E of the log-normal 
distribution. In Cooray & Milosavljevic (2005b), we found 
E ~ 0.23 to describe the field-galaxy luminosity function 
in the K-band (Huang et al. 2003). Here, to describe 2dF- 
GRS data in the fej-band, we found a lower dispersion, with 
E ~ 0.17. While the exact reason for difi'erences between 
the dispersion at two wavelengths is not understood, since 
E reflects the scatter in the Lc{M) relation, we expect this 
relation for luminosities measured in the bj band to have 
less scatter than in the K-band. A comparison of Yang et al. 



(2005) data shown in Fig. 2 and Lin et al. (2004) data shown 
in Fig. 1(a) of Cooray & Milosavljevic (2005b) suggests this 
may be the case, but reasons for the difference in scatter is 
yet to be understood. Incidently, a value for the dispersion 
of 0.17 is in good agreement with the value of 0.168 found 
for the dispersion of central galaxy luminosities by Yang et 
al. (2003b), where these authors used a completely different 
parameterization for the CLF then the one described here. 

In Eq.0 the normalization factor fc{M) in the central 
galaxy CLF captures the efficiency for galaxy formation as a 
function of the halo. In Cooray & Milosavljevic (2005), this 
was set by requiring J ^{L\M)LdL equals the average total 
luminosity of galaxies Ltot(M) in a halo of mass M. This 
condition does not include the fact that at low mass halos, 
galaxy formation is inefficient and not all dark matter halos 
host a galaxy. This is equivalent to modifying the halo mass 
function at the low-mass end to select only halos that host 
a galaxy. Motivated by the halo occupation number models 
for central galaxies (e.g., Kravtsov et al. 2004), where not all 
low mass halos occupy galaxies, to fit the low-end luminosity 
data of the 2dFGRS galaxy LF, we allow a description of the 
form 



MM) 



1 + eii 



f log(M) - log(Afmin 



(2) 



with parameters A/min = 5 x 10^" Msun, and cr = 0.75. These 
parameters were determined by comparing the 2dFGRS LF 
at the low-end as discussed in Section 3. This efficiency func- 
tion is such that it is 0.1 when M ~ 10^*^ but is unity 
when M > few times 10^^ Mq. When describing the envi- 
ronmental LFs, we will continue to use this form. 



2.2 Satellites 

For satellites, the normalization A{M) of the satellite CLF 
can be obtained by defining Ls{M) = Ltot(M) — Lc{M) and 
requiring that Ls(M) = /^'"'" $s{L\M)LdL with /a(L) = 1, 
where the minimum luminosity of a satellite is Lmin- In the 
luminosity ranges of interest, our CLFs are mostly inde- 
pendent of the exact value assumed for Lmin, as long as it 
lies in the range (10^ — W^)Lq. For the maximum lumi- 
nosity of satellites, following the result found in Cooray & 
Milosavljevic (2005b), by comparing predictions to the K- 
band cluster LF of Lin & Mohr (2004), we set Lmax = Lc/2. 
A comparison to 2dFGRS CLFs as measured by Yang et al. 
(2005), however, suggested that such a sharp cut-off is in- 
consistent and that to account for scatter in the total galaxy 
luminosity, as a function of the halo mass, one must allow 
for a distribution in Lmax- Instead of additional numerical 
integrals, we allow for a luminosity dependence with the 
introduction of fs{L) centered around the maximum lumi- 
nosity of satellites such that $s(L|Af) does not go to zero 
rapidly at Lmax. By a comparison to the data, we again 
found a log-normal description with 



ML) = - 



1 + erf 



log(Lc/2.0) 



(3) 



where as = 0.3. The description here is such that /l = 1 
when L < Lmac = Lc/2, but falls to zero at a luminos- 
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Figure 3. The conditional luminosity functions <1?(L|A/) for a variety of mass bins as labeled in each of the panels (with ranges written 
in logarithmic values of the halo mass). The data are measurements by Yang et al. (2005) using the galaxy group catalog of the 2dFGRS 
survey. To describe conditional luminosity functions, 7(A/) varies from —1 at the high mass end to at the low mass end. The comparison 
suggests the evidence of log-normal description for central galaxies and provides an overall better fit to the data than the description 
considered by Yang et al. (2005) based on a priori assumed Schechter (1976) functions for the CLF. Our model, however, under predicts 
the luminosity of central galaxies in halos of mass ~ 10^2 M0. While paramet ers in our model can be modified to fit data better, especially 
a modification to the Vale & Ostriker (2004) Lc{M) relation, we have not pursued this possibility as the overall description is adequate 
enough for the statistical study considered here. Furthermore, it is also unclear to what extent halo mass estimates at the low-mass end 
of the 2dFGRS galaxy group catalog can be trusted and the extent to which selection biases contaminate the CLF determinations at 
this low mass end. 



ity beyond Lc/2 avoiding the sharp drop-off at Lc/2 with 
/s(L) — 1. When model fitting the LF, or the LF as a 
function of the environment, this description is unimpor- 
tant as the central galaxies dominate galaxy statistics. This 
is due to the fact that, as discussed in Cooray & Milosavl- 
jevic (2005b), the LF is dominated by central galaxies in- 
stead of satellites, which has also been noted by Zheng et al. 
(2004) when describing the stellar baryonic mass function. 
Though fs{L) does not matter for the LF, it is important 
when comparing the CLF of galaxies in halos with a nar- 
row mass range and when the CLF is measured with central 



galaxies removed from the data. With regards to satellites, 
note that Ltot(M) equals Lc(M) for M < 10"Mq, and thus 
Ls{M) = O at low halo masses. Such halos only have a single 
galaxy, at the center of the halo. 

In Eq. Q j{M) is taken to be a function of the halo 
mass. In Cooray & Milosavljevic (2005b), we found 7 ~ — 1 
based on model fits to the cluster LF of Lin et al. (2004), 
and we used that value in all mass scales there when mod- 
eling the K-band LF. Here, based on a comparison to CLFs 
measured in the 2dFGRS galaxy group catalog, we find that 
7(M) is, in fact, a function of mass that varies from -1 at 
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cluster scales to at masses corresponding to groups with 
few galaxies. A similar mass dependent variations on the 
faint-end slope, parameterized by a{M) in Schechter func- 
tion forms for the CLF was found by van den Bosch et 
al. (2005) when model fitting the CLF using 2dF LF and 
luminosity-dependent galaxy bias. While the differences as- 
sociated with variations to the slope, as a function of mass, 
are significant when CLFs between galaxy clusters and poor 
groups are considered, in the case of the LF, which averages 
over CLFs of various mass scales, the difference resulting 
from either assuming a mgiss-dependent slope for 7 or an 
average value, such as 7 = —0.5, is minor. Again, this is 
due to the fact that the LF is dominated by central galaxies 
rather than the satellites. Thus, we will set 7(A'/) ~ —0.5 
when calculating the LF and the environmental LFs in the 
present paper, but when describing the CLFs of Yang et al. 
(2005), we allow for a mass dependence for 7 (see, Figure 3). 



2.3 Galaxy types 

Here, in addition to the 2dFGRS LF, we will also model 
the LF of galaxy types, broadly divided in to two classes 
involving red, or early-type, and blue, or late-type, galaxies. 
A previous modeling of galaxy types using CLFs is described 
in van den Bosch et al. (2003); our approach difi^ers because 
of the overall division of the sample to central and satellite 
galaxies. Thus, the division to two types applies to these two 
components, separately, in a given halo. As it is clear, the 
CLFs we have developed facilitate this easy separation. 

As we find later by comparing to 2dFGRS data, central 
galaxies tend to be early- type when found in massive halos, 
corresponding to groups and clusters, but late-type when 
in low mass halos with a few or no satellites. Thus, the 
division of central galaxies to the two types can simply be 
described as a function of mass. While we have extracted this 
division in a function form here, we have not investigated the 
underlying reasons how galaxy types in the center of halos 
change from primarily late to early type, as the halo mass 
function increases. The form is given and it remains to be 
seen if this analytical description is reproduced in rmmcrical 
or semi-analytic models of galaxy formation or not. 

In the case of satellites, a comparison to the 2dFGRS 
CLFs measured by Yang et al. (2005) shows that the divi- 
sion to early- and late- type galaxies is simply not a function 
of mass alone, but rather a function of both mass and lumi- 
nosity of the galaxy. For example, in low mass halos corre- 
sponding to galaxy groups, the low luminous satellites with 
luminosities less than 10'^ Lq tend to be mostly late type or 
blue galaxies, while the bright satellites, with luminosities 
around 10^" Lq are dominated by early type or red galaxies. 

To describe these general behaviors, we consider the 
division of central and satellite conditional luminosity func- 
tions to early- and late-type separately, and write 



*early-cen(I'|M) 
*early-sat(i|M) 



-l'c(I/|M)/early-cen(M) 
*3(L|M)/early-«at(M,L), 



(4) 



where the two functions that divide between early- and late- 
types are taken to be functions of mass, in the case of cen- 
tral galaxies, and both mass and luminosity in the case of 



satellites. Note that these fractions arc defined with respect 
to the total galaxy number of a halo. Since the early- and 
late-type fractions should sum to unity, late-type fractions 
are simply [1 - /oariy-ccn(M)] and [1 - /eariy-sat (M, L)] for 
central and satellite galaxies, respectively. Given this simple 
dependence, we do not write late-type fractions separately. 

In the case of central galaxies, we assume that there is a 
smooth transition between a dominant fraction of late-type 
to a dominant early-type fraction, as a function of the halo 
mass, as the halo mass is increased. A description that fits 
the 2dFGRS CLFs of Yang et al. (2005) was 



/early— cen(-^) — 



1-l-erf 



log(M) - log(Mc, 

(Nearly— cen 



, (5) 



with Mceti = 5 X 10 Mq and fTeariy-cen = 2.0; As the 
halo mass increased, a large width for the transition from 
predominantly late-type to early-type galaxies in the centers 
of massive halos results in an early-type fraction that never 
falls to unity or zero either of the two ends in the mass 
ranges of interest. 

For satellites, at high mass halos, we found early type 
fraction to be roughly two-thirds while at the low mass end, 
this fraction decreases to around one-third. At the low end 
of halo masses where satellites are found, however, this frac- 
tion is luminosity-dependent. The model that describes this 
behavior is 



/early-sat (M,L) = ^ff(M) + ^h{L) + ^ , 

where. 



(6) 



9iM) 
h{L) 



1 + erf 
1 + evi 



f log(M) - log(M3at) 

\ CTaat 

\0g{L) - l0g(Lsa 
Csat 



(7) 



with Msat = 3x 10^^ M0, L^at = 3x 10^ L^un, and cTsat = 1- 
This function varies between 1/3 and 2/3 when low- mass, 
low-luminosities to high-mass, high-luminosities. Note that 
satellites are subhalos that have merged with a central halo. 
Thus, before these galaxies became satellites, they were in 
fact central galaxies. Thus, the fact that the fraction of early- 
to-late satellite galaxies is both mass and luminosity depen- 
dent should not be considered a drawback in this description 
or the halo approach to galaxy statistics in general. In fact, a 
model where redshift dependences are also included, includ- 
ing a merging hierarchy, one may able to start with simple a 
description for early-to-late type galaxy fraction that is mass 
dependent alone and understand how the merging of these 
galaxies result in the fractional dependence of early galax- 
ies, say, as satellites in massive dark matter halos. Clearly, 
such work follows underlying motivations of semi-analytical 
models of galaxy formation. Here, we provide the relations 
that needed to be explored in such an approach. 

In Figure 1, as a summary, we show, as a function of the 
halo mass various fractions encountered in our model. Note 
that the central galaxy fraction falls below unity at halo 
masses below 10^^ Mq, which is due to the efficiency factor 
we included to account for the fact that not all halos at the 
low-end may host a galaxy. Setting /c(M) = 1 results in an 
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Figure 4. The galaxy LF in tlie 2dFGRS survey. Tlie measured data points are from Croton et al. (2004). In addition to the total 
luminosity function, we also show the early (red lines) and late-type (blue lines) LPs. The dotted lines are the central galaxies while 
the dashed lines are for satellites. In (a), we set fc{M) = 1. This results in a faint-end slope of — ~ 1.3, which is steeper than the 
measured value of ~ —1.05. In order to account for the flattening, in (b), we allow for /c(M) to be mass dependent with the form given 
in equation 121 



over prediction for the abundance of galaxies at that low- 
luminosity end. We discuss this in the context of modeling 
the 2dFGRS galaxy LF (Section 3). 



2.4 Central and total-luminosity relation 

The main ingredient in the modeling the LF using this 
empirical approach is the Lc(M) relation. As discussed in 
Cooray & Milosavljevic (2005b), the shape of the Lc{M) 
relation determine the shape of the LF; The slope of this 
relation is directly reflected in the faint-end slope of the LF, 
while the scatter of this relation determines the exponential- 
like drop off of the LF at the bright-end. 

For Lc{M) relation, here we make use of the suggested 
relation in Vale & Ostriker (2004) . These authors established 
this relation by inverting the 2dFGRS luminosity function 
given a analytical description for the sub-halo mass func- 
tion of the Universe (e.g., De Lucia et al. 2004; Oguri & 
Lee 2004). The relation is described with a general fitting 
formula given by 



L(M) = Lo 



{M/Mi 



(8) 



For central galaxy luminosities, the parameters are Lq = 
5.7x 10® Lq, Ml = 10"Mq, a = 4.0, b = 0.57, c = 3.72, and 
d = 0.23 (Vale & Ostriker 2004). These values are different 
from Cooray & Milosavljevic (2005b) since we use K-band 
luminosities there, while the relation given in this paper is 
expected to describe the 2dFGRS data adequately, as it is 
extracted from 2dFGRS LF in the 6,7-band of Norberg et al. 
(2002a). 

For the total galaxy luminosity, as a function of the 
halo mass, we also use the fitting formula in equation 0, 
but with c — 3.57, which we picked based on model fits to 



the 2dFGRS CLFs of Yang et al. (2005). At the massive 
end, the total luminosity can alternatively be described by 
a power-law or a double power-law with the break around 
~ 3.5 X 10" Mq (Vale & Ostriker 2004). The constructed 
LF does not change if power-law behavior is enforced at the 
high-end of the halo masses. This is because the average LF 
is dominated by central galaxies on any scale. The overall 
shape of the LF is strongly sensitive to the shape of the Lc- 
M relation, and it's scatter, and less on details related to 

the Ltot M relation. 

In Figure 2, we show the two relations used for central 
galaxy and total galaxy luminosity of a given halo, on av- 
erage, as a function of the halo mass. Incidently, we found 
the relation extracted by Vale & Ostriker (2004) to be in 
good agreement with direct measurements of central galaxy 
luminosities in the 2dF galaxy group catalog by Yang et al. 

(2005). Similarly, the Ltot M relation based on Vale & 

Ostriker (2004) agrees with total luminosity measurements 
independently obtained by a different technique, involving 
clustering relation, in van den Bosch et al. (2005). In fact, 
these agreement are remarkable and suggest that we have a 
good starting point to build up a model for the galaxy LF 
both as a function of the environment and galaxy type. 



2.5 Conditional Luminosity Functions: A 
Comparison to 2dF 

Now that we have an analytical description for the CLF with 
parameters determined either by results already in the lit- 
erature, such as the Lc{M) from Vale & Ostriker (2004), or 
based on model fitting the data, we can discuss how well our 
models fits the 2dF CLFs of Yang et al. (2005). These CLFs 
were previously modeled with a priori assumed Schechter 
(1976) functions following the models in Yang et al. (2003b), 
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but here we make use of the log-normal description for cen- 
tral galaxies and power-laws for satellite galaxies to build 
the CLF. 

In Figure 3, we show a comparison of our CLFs to those 
extracted from 2dF galaxy group catalog by Yang et al. 
(2005). We show our models in mass ranges they measured 
the CLFs. These CLF model fits can in fact be compared 
with Figures 9 and 11 in Yang et al. (2005). The models 
fits generally support our log-normal description for cen- 
tral galaxies and the separate description of power-law for 
satellite galaxies. The model fits require that the slope of the 
power- law be 7 ~ —1 in cluster scales, as found by Cooray & 
Milosavljevic (2005b) by comparing to the cluster luminos- 
ity function of Lin et al. (2004) in the K-band, but flattens 
to 7 ~ at galaxy group scales. The transition from late- 
to-early type is also adequately modeled with our simple 
description, except that we find our models to under pre- 
dict the number of bright, and potentially central, galaxies 
at poor galaxy group mass scales. This difference may be 
modeled by updating the Lc(M) at these mass scales, but 
given the overall adequate description, and the fact that we 
wanted to build this model with the least number of param- 
eter variations as possible but by using existing results from 
the literature (such as from Vale & Ostriker 2004), we have 
not pursued such a possibility here. It is also not clear how 
well masses have been estimated by Yang et al. (2005) for 
low mass galaxy groups where few galaxies are found in each 
group. 

As discussed in Yang et al. (2005) for 2dF b-band data, 
and in Cooray & Milosavljevic (2005b) for the K-band data, 
the CLF represents galaxy statistics better than the LF 
when wide-field data sets are available in which redshifts 
are measured for tens of thousands of galaxies. While Yang 
et al. (2003b) described CLFs using a priori assumption that 
<I>(L|M) is given by the Schechter (1976) function with no 
separation to central and satellite galaxies, as it is clear from 
Figure 3, our description involving log-normal distribution 
for central galaxies and a power-law for satellite galaxies may 
provide a better description. However, the peak of galax- 
ies at the bright-end may be due to a problem associated 
with mass assignment in the 2dFGRS galaxy group cata- 
log by Yang et al. (2005; van den Bosch, private communi- 
cation). The same reasons could also explain the apparent 
increase in bright galaxies at low mass halos, such as poor 
groups, when compared to our model predictions. Even if 
improvements are minor, when compared CLFs based on 
the Schechter form used in Yang et al. (2005), the method 
suggested here based on a division of the galaxy sample to 
central and satellite galaxies is more physical. As described, 
the model considered also provides a more useful approach 
to divide galaxies to galaxy types, which is advantageous 
since we are trying to get an understanding of how galaxy 
types are distributed in varying dark matter halo masses. 



3 LUMINOSITY FUNCTIONS 

Given our model for the CLFs, we can now construct the 
LF, which is an average of CLFs in mass with the halo mass 
distribution given by the mass function. Here, we use the 




Mass (h-^ M^J 

Figure 7. The probability distribution function of halo masses 
hosting 2dF galaxies, when averaged over the LF. The black lines 
are the total galaxy sample, while red and blue lines show the 
sample divided to early- and late-type galaxies following mea- 
surements by Croton et al. (2004). 

Sheth & Tormen (1999; ST) mass function dn/dM for dark 
matter halos. This mass function is in better agreement with 
numerical simulations (Jenkins et al. 2001), when compared 
to the more familiar Press- Schechter (PS; Press & Schechter 
1974) mass function. 

Given the mass function, the galaxy LF is 

^^{L) = ML\M)^dM , (9) 

where i is an index for early and late type galaxies. The 
conditional luminosity function for each type involves the 
sum of central and satellites. To understand our model for 
the LF, we will plot these two divisions, as well as the sum, 
separately in each of the LF figures shown in this paper. 

Figure 4 presents a general fit to the 2dF LF, where in 
Fig. 4(a), we assume /c(M) = 1. While the model describes 
data adequately at magnitudes below -19, the fit at the low- 
end of luminosity is poor. Our model suggests that the low- 
end slope of the LF should be around —1.33 to —1.25, with 
the latter coming from our assumption that Lc{M) ~ 
at low mass halos (Cooray & Milosavljevic 2005b). The fiat- 
tening of the slope at the low-end of luminosity, as measured 
in the 2dF LF, may be a refiection of either (1) galaxy se- 
lection function such that faint galaxies are missed, or (2) 
a real effect in the Universe such that low mass halos do 
not host a large number of central galaxies. The difference 
from the model in Fig. 4(a), when compared to data, can 
be reduced with the inclusion of a mass-dependent /c(M) 
function. In Fig. 4(b), we show the case with equation |21 
Note that the underlying reasons for this mass function can 
either be a selection effect or a real effect. We cannot dis- 
tinguish between the two possibilities; but, it is likely that 
most selection biases are already accounted when construct- 
ing the LF and the effect we are seeing mostly is due to the 
fact that not all low mass dark matter halos host a galaxy. 
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Figure 5. The mass dependence of the 2dF LF. The plotted data and hne styles are same as in Figure 4. Here, from (a) to (d), we show 
the contribution to the luminosity function from a galaxies in the mass ranges 10^ Mq to 10^^ ^0' lO'^^ Mq to lO'^^ ^0' 10^^ Mq 
to 10^'* Mq, and 10^'* Mq to 10^® Mq respectively. Central galaxies in low mass halos determine the shape of the LF at the faint-end, 
while bright galaxies in massive halos, such as groups and clusters, determined the LF at the bright-end, where an exponential-like cut 
off is observed. 



To understand the mass dependence of the LF, in Fig- 
ure 5, we plot the LF separated in to mass bins between 10^ 
Mq to 10^^ Mq. Galaxies in low mass halos dominate the 
statistics of the LF at the faint-end; in fact, 2dF luminosity 
function at magnitudes fainter than -19 is associated with 
galaxies in halos with masses below 10^^ Mq. On the other 
hand, the exponential decrease in the Schechter form for the 
LF at the bright end is associated with galaxies in halos of 
meiss lO'^^ Mq and above. As shown in Figure 1, the Lc(M) 
relation begins to flatten at halo masses above lO" Mq. 
The characteristic luminosity, , can be identified with the 
luminosity of galaxies in halos of this mass scale. The ex- 
ponential drop-off, instead of a sharp-cut off is a reflection 
of the scatter in the L(M) relation above this mass scale as 
was explained in Cooray & Milosavljevic (2005b). To fit the 
2dF LF, we require E ~ 0.17; this is lower than the value 



of ~ 0.23 found in Cooray & Milosavljevic (2005b) when 
modeling the K-band LF of Huang et al. (2003). 

To understand the relative distribution of mass given 
the galaxy luminosity, we also calculate the conditional 
probability distribution P{M\L) that a galaxy of given lu- 
minosity L is in a halo of mass M as (Yang et al. 2003b): 



P{M\L) dM = 



^{L\M) dn 
$(L) dM' 



dM . 



(10) 



In Figure 6, we show these conditional probability distribu- 
tion functions, as a function of the halo mass, when L = 10* 
h~^ Lq to lO^'^ Lq. These probabilities show a peak 
at low masses, associated with central galaxies, and a tail 
towards higher masses, associated with satellites of the same 
luminosity. As the mass scale is increased, the peak related 
to central galaxies broaden since the Lc(M) relation in- 
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Figure 6. The conditional probability distribution function of halo mass P(M\V) to host a galaxy of the given luminosity, as function 
of the halo mass. The black lines are the total galaxy sample, while red and blue lines show the sample divided to early- and late-type 
galaxies. These are the same models that were used to describe the LFs of Croton et al. (2004) and shown in Figure 4. Prom (a) to 
(d), we show these probabilities for L = 10*/i-2 Lq, L = lO^h-^ Lq, L = lO^^h'^ Lq, and L = lO^^h''^ Lq, respectively. The low 
mass end peak, which tends to be narrow for lower luminosities, are the central galaxies, while the tail to high masses, is associated with 
satellite galaxies. The width of the central galaxy luminosity peak increases since the Lc{M) relation increases while flattening such that 
one encounters a fractionally higher mass range with an increase in luminosity. Low luminous late-type galaxies are found in low-mass 
halos, while luminous red galaxies are in halos with masses corresponding to galaxy groups and clusters. 



creases slowly with increase in mass such that one encoun- 
ters a fractionally higher mass range with an increase in 
luminosity. Fainter late-type galaxies are in low mass halos, 
while brighter early type galaxies are in galaxy groups and 
clusters. 

By integrating the conditional probability distribution 
functions over luminosities, we can calculate the the proba- 
bility distribution of mass associated with the 2dF LP (van 
den Bosch et al. 2003): 



P{M) dM = 



/ ^{L\M)dL dn 
J ^{L)dL 



dM 



dM . 



(11) 



As shown in Figure 7, the LF is dominated by galaxies that 
occupy dark matter halos in the mass ranges between lO'^" 
Mq to 10^^ While central galaxies dominate statistics 



in this mass range, satellite galaxies become the dominant 
contributor to the LF from each mass scale. This, however, 
does not imply that at each luminosity, satellites dominate, 
but rather, as Figure 4 shows, central galaxies dominate. 
Note that the behavior of central galaxies is similar to the 
average, but on the other hand, early type or red galaxies, 

l0. 



are primarily in halos with masses above 10^^ Mr; 



3.1 Galaxy Bias 

Another useful quantity to compare with observed data is 
the galaxy bias, as a function of the luminosity. Using the 
conditional LFs, we can calculate these as 
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Figure 8. The environmental LFs, or ^(iliSgai), wliere (5gai is the galaxy overdensity measured over a volume of 8 h''^ Mpc. The 
measured data are from Croton et al. (2004). From (a) to (f), the overdensity varies from ~ —0.9 to greater than 6. In each of the panels, 
we show the "best-fit" models for the conditional mass function, dn{M\5^^i) / dM , based on parameters A{&g^i) and «((5gai) following 
Eg. 1141 The line styles are same as Figure 4. For comparison, we also show the average LF in Figure 4(b), with dot-dashed lines. 
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Figure 9. (a) The predicted mass functions, as a function of the environment captured by the galaxy overdensity over a size scale of 8 
Mpc. The dot-dashed line is the standard ST mass function, while other curves show the mass function from voids with (5gai < 
to dense regions with iS > 1. (b) The ratio of environmental mass functions to the mean mass function. Note the sharp cut-off at high 
masses when S < —0.5. In the case of dense regions when 5 2> 1, the mass function increase at all mass scales, relative to the mean, but 
with a greater increase at high masses. 



(12, 

where &haio(Af) is the halo bias with respect to the linear 
density field (Sheth, Mo & Tormen 2001; also, Mo et al. 
1997) and i denotes the galaxy type. In Figure 11, we show 
the galaxy bias as a function of the luminosity. We also di- 
vide the sample to galaxy types. 

While the bias factors have similar shapes, early-type 
galaxies are biased higher relative to the late-type galaxies; 
the difference in the bias factor between early- and late-type 
galaxies is at the level of ~ 5%. At low luminosities, the total 
sample bias is close to that of late- type galaxies, while at the 
bright-end average bias for the whole sample is close to that 
of early-type galaxies. We also note an important differences 
in the bias when comparing satellite galaxies and central 
galaxies; Satellite galaxies, on avearge, have a higher bias 
factor at all luminosities when compared to central galaxies. 
This is due to the fact that satellite galaxies are preferen- 
tially in higher mass halos that are, on averge, biased higher 
with respect to the linear density field. The average bias 
factor, for the whole sample, however is dominated by cen- 
tral galaxies, for same reasons the LF is also dominated by 
central galaxies. Bias measurements as a function of galaxy 
type exist in the form of clustering information such as the 
correlation length as a function of luminosity and type (Nor- 
berg et al. 2002b). While we have not attempted to convert 
this clustering information to obtain bias fe(L) as a function 
of galaxy type here, since it involves an additional step of 
modeling, our next improvement in this approach is to com- 
pute clustering statistics, in which case a direct comparison 
could easily be made. This clustering information has been 
used in van den Bosch et al. (2003) when modeling the CLFs 
appropriate for the 2dFGRS survey. 




Figure 11. Galaxy bias as a function of luminosity calculated 
from conditional luminosity functions {solid line); we include the 
separate contributions from central galaxies (dotted line) and 
satellites (dashed line). Also shown are the SDSS of Zehavi et al. 
(2004), and a fit to the galaxy bias in 2dF data of Norberg et al. 
(2002b) (dot- dashed line). We also show the bias for galaxy types 
(early and late). Late- type galaxies are expected to be in low- 
dense regions and their bias factor, relative to early type galaxies, 
would be lower. Satellite galaxies, regardless of the type, are in 
more massive halos and, thus, have higher bias factors relative to 
central galaxies. The average bias factor, for the whole sample, 
however is dominated by central galaxies, for same reasons the 
LF is also dominated by central galaxies. 
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Figure 10. The probability distribution function of halo masses hosting 2dF galaxies as a function of the environment as measured by 
the galaxy overdensity in 8 Mpc volumes. Here, we plot P(log M|(5gai). The black lines are the total galaxy sample, while red and 
blue lines show the sample divided to early- and late-type galaxies following measurements by Croton et al. (2004). Form (a) to (d), 
we show four ranges in corresponding to both under dense and over dense regions. In the case of under dense regions, halos with 
mass 10^^ Mq host most galaxies and the the fraction of galaxies is dominated by the late types (blue galaxies). In the case of over 
dense regions, a variety of mass scales contribute to the luminosity function of that region, with early-types dominating the fraction from 
massive halos, while late-types dominate the fractional contribution from low mass halos. 



4 ENVIRONMENTAL LUMINOSITY 
FUNCTIONS 

Having discussed average statistics, we now focus on the 
LFs measured by Croton et al. (2004) as a function of the 
galaxy overdensity, (5gai. These overdensities correspond to 
a volume of radius 8 Mpc, and are measured following 

5gal = {Ng - Ng)/Ng SUch that -1 > (5gal. 

To calculate the environmental LFs given by $(L|5gai), 
the LF given 5gai, following Mo et al. (2004), we make one 
important assumption; We assume that the CLFs, $(LjAf) 
are independent of the environment. The motivation for such 
an assumption is inherent in the halo model (Cooray & 
Sheth 2002), where one makes the assumption that galaxy 
distribution can be associated with the halo mass rather 



than the environment. Similarly, observational studies based 
on the Sloan survey indicate that galaxy color may also be 
independent of the environment, implying that the CLFs for 
galaxy types are not dependent on 5gai (Blanton et al. 2004). 
Thus, given that "I>(L|M) is (5gai independent, all variations 
in <E>(L|5gai) must come from variations in the halo mass 
function, as a function of 5gai. The halo mass function, in 
fact, is (5dm dependent; the peak background split (Sheth & 
Tormen 1999,2002) makes use of this dependence to extract 
information on, for example, relative biasing of the halos, as 
a function of mass, to calculate halo bias factors. 

Thus, to calculate the LF as a function of the environ- 
ment, as defined by the galaxy overdensity, we define 

ML\S,.i) = J^^ ML\M)^^^^^^dM , (13) 



© 0000 RAS, MNRAS 000, 000-000 



14 Cooray 



where dn{M\5gs_i) /dM is the conditional halo mass func- 
tion. Unfortunately, analytical techniques are not adequate 
enough to reliably calculate the conditional mass function, 
fin/drfi(Af|5gai) given the galaxy overdensity, 5gai (see, e.g.. 
Mo ct al. 2004 and discussion below). Here, making use of 
the assumption that CLFs are independent of iJgai, we use 
observed measurements of $(L|(5gai) to extract information 
on dn{M\5ga,i)/dM, the conditional mass function of daxk 
matter halos as a function of the galaxy overdensity. In this 
approach, dn{AI\5^ai) / dA4 necessary to reproduce Croton et 
al. (2004) data can potentially be compared with either im- 
proved models of the conditional mass function. While this 
comparison is beyond the scope of this paper, as it involves 
understanding certain aspects of the mass function, we plan 
to return to this topic later. Later in this Section, we will, 
however, comment on the void mass function or the condi- 



tional mass function when Seed 



-1. Since Croton et al. 



(2004) measurements include the LF of galaxies in voids, 
we can directly establish the void mass function using our 
modeling technique. 

To extract information on dn(M|5gai)/dM, motivated 
by numerical simulation based measurements in Mo et al. 
(2004) of the conditional mass function, we assume one can 
model this as 



dn{M\5,^i)/dM = A(5) 



1 + 



M 

1013 



dn/dM 



(14) 



where A{5) and a{5) are two parameters we will extract from 
the data given the average mass function dn/dM following 
the ST mass function. 

These models fits to the Croton et al. (2004) data are 
shown in Figure 8. For comparison, we also show the av- 
erage luminosity function (shown in Figure 4). Our model 
fits generally describe the data, except in overdense regions, 
such as when (5gai > 6, when we underpredict the abun- 
dance of early type galaxies and overpredict the abundance 
of late-type galaxies. 

The conditional mass functions that describe these envi- 
ronmental luminosity functions are shown in Figure 9, where 
we plot separately the mass functions and the ratio of mass 
functions to the average ST mass function. Under dense re- 
gions show a sharp decrease in the abundance of high mass 
halos, while the over dense regions are such that it is a scaled 
version of the ST mass function, but with an increasing 
abundance of high mass halos. In principle, these mass func- 
tions should be described by the conditional mass functions 
that are used to calculate merger trees of dark matter ha- 
los or based on the peak-background split (Shcth & Tormcn 
1999, 2002). However, the require modification to the mass 
function with S ^ Sc- 5L,m and cr^(M) a'^{M) - a^{M') 
did not produce shapes of the extracted mass functions. In 
this replacement, Sc ~ 1.686 is the standard overdensity for 
collapse, M' — A/'iirR"^ pm{l + dm), and 5m is the overden- 
sity in mass that corresponds to the overdensity in galaxies, 
and SL,m is the linearized overdensity corresponding to the 
mass overdensity. The main reason for the difficulty in ob- 
taining dn{5gsi)/dm is that we do not yet have an accurate 
model for the relation between (5gai, which is non-linear, and 
5m and 5L,m- In fact, these relations capture the biasing of 
galaxies with respect to the linear density field. In Mo et 



al. (2004), authors used simulations to estimate the envi- 
ronmental mass functions, as a function of the galaxy over- 
density, which were then compared with same LFs directly 
measured in mock catalogs. Here, we use observed data to 
extract information on the conditional mass functions. We 
find reasonable agreement with the mass functions plotted 
in their Figure 2 ^ 

To further understand the dependence of these LFs on 
halo masses, we calculate conditional probability for halos 
of mass M to host galaxies given the environmental over- 
density, (5gai. These probabilities arc calculated from 



/ ^L\M)dL dn(M|^gai) 
/ ^L)dL dM 



(15) 



In Figure 10, we plot P(log M|5gai). These probabilities show 
mass scales that arc important for the environmental LF, 
as a function of the galaxy overdensity. In addition to the 
total distribution, we also show probabilities in terms of 
the galaxy type, and separated to both central and satel- 
lite galaxies. The halos that contribute to the LF in under 
dense regions host a higher fraction of late-type galaxies rel- 
ative to early type ones. On the other hand, in over dense 
regions, one finds both late-type and early-type galaxies, 
with a large fraction of early-type galaxies coming from more 
massive halos or galaxy groups and clusters, while the late- 
type fractional contribution is dominated by galaxies in low 
mass halos. 

To compare with measurements in Croton et al. (2004), 

wo calculate one more quantity involving the moan luminos- 
ity per galaxy, as a function of the galaxy density environ- 
ment: 



r ^L\5g,i)LdL 
Ng~ , $(L|5gai)dL 



PL 



(16) 



Here, I/min was set at an absolute magnitude of —17 follow- 
ing the measurements of Croton ct al. (2004). The Croton 
et al. (2004) data and the same quantity based on model 
fits in Figure 8 are plotted in Figure 12. They show general 
agreement, except in high dense environments, the mean 
luminosity per galaxy, in the case of early type galaxies, 
is some what higher in our models when compared to the 
data. This is due to the fact that our models under predict 
the abundance of early type galaxies in dense environments, 
making them brighter on average than observed. 



4.1 Galaxies in Voids 

As a final use of data in Croton et al. (2004), we use their 
environmental LF in extreme voids, with — 1 < 5gai < —0.9, 
to extract information on the conditional mass function of 



1 This is true only if the labels in Figure 2 of Mo et al. (2004) 
is reversed from what is labeled there. As shown in our figure 10, 
under dense mass functions show a low abundance of halos rather 
than the increase abundance, relative to the over dense regions, 
as suggested in Figure 2 of Mo et al. (2004). This is likely to be 
a misprint. 
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Figure 12. The mean luminosity per galaxy as a function of the 
environment, as measured by (5gai. The data plotted are direct 
measurements established by Croton et al. (2004) based on the 
measured environmental LFs, while the lines show the predictions 
based on our model fits shown to describe the same environmental 
LFs (Figure 8). The solid circles show the total galaxy sample, 
red squares show early type galaxies, and blue triangles show late- 
type galaxies. We find a general agreement, except at the high 
dense regions, where the mean luminosity per galaxy for early 
type galaxies is bit higher in our models when compared to the 
data. This is due to the fact that in our models, we underpredict 
the abundance of early type galaxies, relative to observed LF. 



the same environment. Our results our summarized in Fig- 
ure 13. Note that the void LF show a sharp decrease in 
the abundance of galaxies with absolute magnitudes brighter 
than -19. Since our conditional LFs are independent of S^^i, 
this puts a strong constraint on the void mass function at 
the high mass end. Since luminosities grow with mass, fol- 
lowing the Lc{M) relation, to restrict galaxy luminosities 
to be fainter than the observed cutoff requires that no ha- 
los with masses greater than 10^'^ Msun be present in voids 
with a fractional abundance Mdn/dM greater than 10~* h'^ 
Mpc"^. 

Given that the void mass function has received some 
attention in the literature, we can make a direct compar- 
ison with our estimate of the mass function with previous 
estimates to the extent analytical formulae are available. 
In Figure 13, we show several comparison drawn from pub- 
hshed attempts to analytically model the void mass function 
(also, see Gottlober et al. 2003). The methods by Goldberg 
et al. (2004) and Patiri et al. (2005) generally fail to describe 
the void mass function both at the low- and high-end of the 
mass function, though the abundance at a halo mass of 10^^ 
Mq is generally produced. The Sheth & van de Weygaert 
(2003) analytic description for the void mass function comes 
closest to describing the mass function required by the 2dF 
extreme void LF. In fact, Sheth & van de Weygaert (2003) 
description managed to capture the low-end turn over in the 
void mass function, which is simply associated with our ef- 
ficiency function fc{M). This model, however, overproduces 



the abundance of halos at the high mass end with a pre- 
diction for the presence of ~ lO'^* Mq halos with the same 
abundance as halos of ~ 10^'^ Mq required by the 2dF LF. 
The dashed lines in Figure 13(a) show the predicted LF 
based on the Sheth & van de Weygaert (2003) void mass 
function. The bright-end is over populated and the turnoff 
in the LF moves to a higher luminosity than seen in the 
data. Our suggestion for the void mass function, based on 
the Croton et al. (2004) extreme void LF, may provide a 
useful guidance to analytically calculate the void mass func- 
tion. As part of an attempt to understand the conditional 
mass functions, we hope to return to this issue again in the 
future. 



5 DISCUSSION AND SUMMARY 

In this paper, we have made use of the scaling relation be- 
tween the central and the total galaxy luminosity of a dark 
matter halo as a function of the halo mass, to construct 
the CLF of galaxies. The CLF provides a powerful tech- 
nique to understand galaxy properties, especially on the 
spatial distribution of galaxies, in the era of wide-field large 
scale structure surveys where large galaxy samples can easily 
be subdivided with adequate statistics (Yang et al. 2003b). 
The CLF is closely related to the halo occupation number, 
which is the average number of galaxies given the halo mass 
(Cooray & Sheth 2002). The halo occupation number cap- 
tures galaxy statistics that simply treat each galaxy in the 
sample equally and is useful when attempting to understand 
statistics such as galaxy clustering statistics, on average. 
Higher order statistics of the galaxy distribution require de- 
tailed statistics of halo occupation beyond the mean. On the 
other hand, galaxies vary in luminosity and color. Therefore, 
a more useful quantity to consider is the conditional occu- 
pation number, ie., the number of galaxies in a given halo 
mass with given luminosity. These conditional occupation 
numbers are in fact the conditional luminosity functions we 
have described here. 

When compared to previous descriptions of the CLF in 
the literature (e.g., Yang et al. 2003b; Yang et al. 2005), we 
make several improvements by dividing the galaxy sample 
to galaxies that are in centers of dark matter halos (central 
galaxies), and in subhalos of a main halo (satellite galax- 
ies). This division is central to are arguments on how the 
CLF is shaped and assumes that different physics govern the 
evolution of these two components. In our underlying phys- 
ical description, following Cooray & Milosavljevic (2005a), 
central galaxies grow in luminosity through dissipationless 
merging of satellite galaxies. These satellites are both early- 
and late-types, with the fraction of early type galaxies taken 
to be dependent on both the halo mass and the satellite 
galaxy luminosity. Underlying physical reasons for this de- 
pendence is not clear, but could very well associated with a 
tidal stripping effect that removes gas and the stellar con- 
tent of late-type galaxies in more massive halos and con- 
vert these galaxies to early-types. It will be helpful to un- 
derstand if numerical and semi-analytic models predict the 
fractions of late to early type galaxies, including the sugges- 
tion that the fraction in satellites changes at a luminosity 
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Figure 13. A comparison of tlie void LF (a) and the void mass function (b). Here, we show the galaxy LF when —1 < 5gai < -l0.9, 
where S^^i is measured over a radius of 8 Mpc. The data points are measurements from Croton et al. (2004), where only the total LF 
in voids was measured. The curves follow the same line-styles as Figure 8, and with dot-dashed lines showing the mean LF in Figure 4. 
In (b), we show the mass function that fits the Croton et al. (2004) data with a thick solid line. The dot-dashed line is the ST rhass 
function. We also show three predictions from the literature on the void mass function. Both Goldberg et al. (2004) and Patiri et al. 
(2005) fail to describe both the low- and high-mass end shapes of the mass function derived here. Interestingly, the analytic description 
by Sheth & van de Weygaert (2004) fits the low-mass end of void mass function, which is captured by our function fc{M) applied to the 
average mass function, but over predicts the abundance of high mass halos in voids. For example, the void LF calculated with the SHieth 
&i van de Weygaert (2004) mass function for voids is shown in (a) with dashed-lines. While the low-luminosity end agrees with the d]ata, 
the high luminosity end, and the turnover in the luminosity function is higher than suggested by the data. This increase is associiited 

with the high mass allows in the Sheth &: van de Weygaert (2004) void mass function. ] 

I 
1 

direct predictions that can compare with observations,; we 
failed to do this due to lack of information on the conditional 
mass function, or the mass function of dark matter hklos 
given the galaxy overdensity. In Mo et al. (2004) predictions 
were made using the conditional mass function measured 
in numerical simulations. Here, we use Croton et al. (2(il04) 
measurements to establish information on the conditional 
mass functions, from environments such as galaxy voidfe to 
dense regions. With these mass functions, we also estinjate 
statistical quantities as probability distribution function of 
halo mass, as a function of the galaxy overdensity. We find 
that the preferred environment of blue galaxies are under- 
dense regions in low mass halos, while the early- type, ; red 
galaxies are mostly in dense environments and dominated 
by satellites of larger mass halos. The shapes of the imass 
functions we have extracted could eventually be compdred 
to numerical simulations or analytical techniques. ; 

Using Croton et al. (2004) LF for galaxies in voids, we 
also establish the void mass function. We find that the void 
mass function is peaked at halo masses around ~ 10^^ ^0- 
Such a peak in the void mass function is predicted in the 
analytical calculation by Sheth & van de Weygaert (2003). 
We do not, however, find a tail to higher halo masses as 
suggested by the description in Sheth & van de Weygaert 

(2003) . The void mass function must sharply turn over; if 
not, one would predict brighter galaxies in void environ- 
ments than suggested by the measurements in Croton et al. 

(2004) . 



around 5 x 10 Lq from more late-types at lower lumi- 
nosities than this value to more early-types. 

Our simple empirical model describes the conditional 
luminosity functions of galaxy types, as well as the total 
sample, measured with the 2dF galaxy group catalog from 
cluster to group mass scales by Yang et al. (2005). A compar- 
ison of our CLFs to ones in Yang et al. (2005; their Figure 9 
and 11) reveal that our model can account for the peak of the 
CLF at high luminosities; this peak is associated with central 
galaxies, which we have modeled with a log-normal distri- 
bution. Similar peaks are also observed in the conditional 
baryonic mass function of stars in semi-analytical models of 
galaxy formation (Zheng et al. 2004) and in the luminosity 
functions of clusters where all galaxies are included (Tren- 
tham & Tully 2002). The slope of the satellite CLF changes 
from -1 at cluster mass scales to at group mass scales. 
Again this could be a refiection of the merging evolution 
of satellite galaxies as hierarchical merging builds up bigger 
halos. It will be interesting to see if a simple extension of 
the analytic model of Cooray & Milosavljevic (2005a), which 
involves dynamical friction, may explain the change in the 
slope as parent halo mass is increased. 

Given the observational measurements of the LF as a 
function of the environment using 2dF data by Croton et 
al. (2004), as measured by the galaxy overdensity measured 
over volumes corresponding to a radius of 8 h~^ Mpc, we 
extended our CLF model to describe environmental lumi- 
nosity functions. While it would have been useful to have 
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To summarize our paper, main results are: 

(1) The galaxy LF, which is an average of CLFs over 
the halo mass function, is primarily shaped by the Lc{M) 
relation; the faint-end slope of the LF reflects the faint-end 
scaling of the Lc{M) relation while the bright-end turn off in 
the LF, generally described by the exponential cut off in the 
Schechter (1976) fitting form, is determined by the scatter 
in the Lc(M) relation. Understanding the Lc{M) relation 
and it's scatter is central to understanding the galaxy LF. 

(2) The fraction of galaxy types as a function of the 
halo mass (Figure 1). The galaxy types are distributed such 
that one finds essentially all late-types in low mass halos 
and mostly early-types in high mass halos. In the case of 
early-types, the fraction is dominated by satellite galaxies 
rather than central galaxies in halo centers. These fractions 
may capture interesting physics such as tidal stripping that 
happens in dense environments and could also be responsible 
for the fractional change of galaxy types, from dominate 
late-type to early-type, as the luminosity of satellite galaxies 
is increased. 

(3) The mass dependent slope for the satellite CLF, 
where "/{M) ~ — 1 at galaxy cluster mass scales and 7 ~ at 
poor galaxy group mass scales. The slope may be a reflection 
of the hierarchical merging process and predictions on 7(M) 
do not yet exist in the literature. 

(4) The conditional mass functions (Figure 9), or the 
mass function of dark matter halos given the galaxy over- 
density measured over a volume corresponding to a radius 
of 8 h~^ Mpc. We will leave it as a challenge to improve an- 
alytical techniques produce the required mass functions to 
compare with what we suggest is needed to explain Croton et 
al. (2004) measurements. Any disagreements, if understood, 
may suggest that our central assumption that the CLF is 
independent of the galaxy overdensity is incorrect. If that's 
the case, galaxy formation and evolution involves additional 
parameters beyond the halo mass and could question the 
viability of the halo approach to describe galaxy statistics. 

(5) Galaxy bias predictions as a function of luminosity 
given the galaxy type (Figure 11). While our predictions 
agree with the bias predicted for total samples, it will be 
useful to understand clustering bias as a function of galaxy 
color as well. In an upcoming paper, we will extend the CLFs 
developed here to describe clustering properties of galaxies 
and a direct comparison to measurements in Norberg et al. 
(2002b). 

Acknowledgments: 

We thank Milos Milosavljevic for his contributions that 
began this project from an attempt to understand a sim- 
ple observed relation (Cooray & Milosavljevic 2005a) and 
Frank van den Bosch for his detailed comments and sugges- 
tions that convinced the author to use this simple relation to 
model the LF. The author also thanks Darren Croton and 
Xiaohu Yang for providing electronic files of various mea- 
surements related from 2dFGRS that are modeled in the 
paper. The author thanks members of Cosmology and The- 
oretical Astrophysics groups at Caltech and UC Irvine for 
useful discussions. 



REFERENCES 

Balogh, M. L., Pearce, F. R., Bower, R. G., & Kay, S. T. 2001, 
Mon. Not. R. Astron. See, 326, 1228 

Barkana, R. & Loeb, A. 2001, Physics Reports, 349, 125 

Benson, A. J., Pearce, F. R., Frenk, C. S., Baugh, C. M., & Jenk- 
ins, A. 2001, Mon. Not. R. Astron. Soc, 321, L7 

Benson, A. J., Lacey, C. G., Baugh, C. M., Cole, S., & Frenk, 
C. S. 2002, Mon. Not. R. Astron. Soc, 333, 156 

Benson, A. J., Bower, R. G., Frenk, C. S., Lacey, C. G., Baugh, 
C. M., & Cole, S. 2003, Astrophys. J., 599, 38 

Berlind, A. A., Blanton, M. R., Hogg, D. W., et al. 2004, 
astro-ph/0406633 

Binney, J. 2004, Mon. Not. R. Astron. Soc, 347, 1093 

Blanton, M. R. et al. 2003, Astrophys. J., 819 

Blanton, M. R., Eisenstein, D. J., Hogg, D. W., & Zehavi, I. 2004, 
astro-ph/0411037 

Bullock, J. S., Kravtsov, A. V., & Weinberg, D. H. 2000, Astro- 
phys. J., 539, 517 

Cole, S. et al. 2001, Mon. Not. R. Astron. Soc, 326, 255 

CoUess, M. et al. 2001, MNRAS, 328, 1039 

Cooray A., Hu W., 2001, ApJ, 554, 56 

Cooray A., Hu W., & Miralda-Escude, J. 2000, Astrophys. J., 
535, L9 

Coo ray, A ., & Sheth, R. 2002, Phys. Rep., 372, 1 

jastro-ph /0206508 I 
Cooray, A., & Milosavljevic, M. 2005, preprint 

jastro-ph/0 503596 I 
Cooray, A., & Milosavljevic, M. 2005, preprint 

jastro-ph/0 504580 I 
Croton, D. J. et al. 2004 Mon. Not. R. Astron. Soc. in press 

llastro-ph/0407537 1 
Dekel, A. 2004, in Multiwavelength Mapping of Galaxy Forma- 
tion and Evolution, eds. R. Bender & A. Renzini, preprint 

( astro-ph /0401503 I 
Dekel, A., & Birnboim, Y. 2004, preprint I astro-ph/0412300 i 
De Lucia, G., Kauffmann, G., Springel, V., White, S. D. M., Lan- 

zoni, B., Stoehr, F., Tormen, G., & Yoshida, N. 2004, Mon. 

Not. R. Astron. Soc, 348, 333 
Dressier, A. 1980, Astrophys. J., 236, 351 

Dressier, A., Oemler, A. Jr., Couch, W. J., Smail, I., Ellis, R. S., 

Barger, A., Butcher, H., Poggianti, B. M., Sharpies, R. M. 

1997, Astrophys. J., 490, 577 
Drory, N., Bender, R., Feulner, G., Hopp, U., Marston, C, 

Snigula, J., & Hill, G. J. 2003, Astrophys. J., 595, 698 
Goldberg, D. M., Jones, T. D., Hoyle, F., Rojas, R. R., Vogeley, 

M. S., & Blanton, M. R. 2005, Astrophys. J., 621, 643 
Goto, T., Yamauchi, C, Fujita, Y., Okamura, S., Sekiguchi, M., 

Smail, I., Bernardi, M. & Gomez, P. L. 2003, Mon. Not. R. 

Astron. Soc, 346, 601 
Gottlober, S., Lokas, E., Klypin, A. & Hoffman, Y. 2003, Mon. 

Not. R. Astron. Soc, 344, 715 
Huang, J.-S., Glazebrook, K., Cowie, L. L, & Tinnoy, C. 2003, 

Astrophys. J., 584, 203 
Jenkins, A., Frenk, C. S., White, S. D. M., Colberg, J. M., Cole, 

S., Evrard, A. E., Couchman, H. M. P., & Yoshida, N. 2001, 

Mon. Not. R. Astron. Soc, 321, 372 
Kay, S. T., Pearce, F. R., Frenck, C. S., & Jenkins, A. 2002, Mon. 

Not. R. Astron. Soc, 330, 113 
Keres, D., Katz, N., Weinberg, D. H., & Dave, R. 2004, preprint 

|astro-ph/04 07095j 
Kravtsov, A. V., Berlind, A. A., Wechsler, R. H., Klypin, A. A., 

Gottlober, S., AUgood, B., & Primack, J. R. 2004, Astrophys. 

J., 609, 35 

Lin, Y., & Mohr, J. J. 2004, Astrophys. J., 617, 879 



© 0000 RAS, MNRAS 000, 000-000 



18 Cooray 



Lin, Y., Mohr, J. J., & Stanford, A. 2004, Astrophys. J., 610, 745 
Mo, H. J., Jing, Y. P., White, S. D. M. 1997, Mon. Not. R. Astron. 
Soc, 284, 189 

Mo, H. J., Yang, X., van den Bosch, F. C, & Jing, Y. P. 2004, 

Mon. Not. R. Astron. Soc, 349, 205 
Norberg, P., et al. 2002a, Mon. Not. R. Astron. Soc, 336, 907 
Norberg, P., et al. 2002b, Mon. Not. R. Astron. Soc, 332, 827 
Og uri, M. fc Lee, J. 2004, Mon. Not. R. Astron. Soc, 

lastro- ph/0401628| 
Pa tiri, S. G., Betancort-Rijo, J. E., Prada, F. 2004, 

|astro-p h/0407513 
Peacock J. A., Smith, R. E., 2000, MNRAS, 318, 1144 
Press, W. H., & Schechter, P. 1974, Astrophys. J., 187, 425 
Rees, M. J., & Ostriker, J. P. 1977, Mon. Not. R. Astron. Soc, 

179, 451 

Schechter, P. 1976, Astrophys. J., 203, 297 
Scoccimarro R., Sheth R., Hui L., Jain B., 2001, ApJ, 546, 20 
Seljak, U. 2000, Mon. Not. R. Astron. Soc, 318, 203 
Sheth, R. K., & Tormen, G. 1999, Mon. Not. R. Astron. Soc, 
308, 119 

Sheth, R. K., & Tormen, G. 2002, Mon. Not. R. Astron. Soc, 
329, 61 

Sheth, R. K., Mo, H. J., & Tormen, G. 2001, Mon. Not. R. Astron. 
Soc, 323, 1 

Sheth, R. K., & van de Weygaert, R. 2004, Mon. Not. R. Astron. 
Soc, 350, 517 

Spergel, D. N., et al. 2003, Astrophys. J. Supp., 148, 175 
Trentham, N., & TuUy, R. B. 2002, Mon. Not. R. Astron. Soc, 
335, 712 

TuUy, R. B., & Pierce, M. J. 2000, Astrophys. J., 533, 177 
Vale, A., & Ostriker, J. P. 2004, Mon. Not. R. Astron. Soc, 353, 
189 

van den Bosch, F. G., Yang, X., & Mo, H. J. 2003, Mon. Not. R. 

Astron. Soc, 340, 771 
van den Bosch, F. C., Yang, X., Mo, H. J. & Norberg, P. 2005, 

Mon. Not. R. Astron. Soc, 356, 1233 
White, S. D. M., & Rees, M. J. 1978, Mon. Not. R. Astron. Soc, 

183, 341 

Yang, X., Mo, H. J., Kauffmann, G., & Chu, Y. Q. 2003a, Mon. 

Not. R. Astron. Soc, 339, 387 
Yang, X., Mo, H. J., & van den Bosch, F. C. 2003b, Mon. Not. 

R. Astron. Soc, 339, 1057 
Yang, X., Mo, H. J., Jing, Y. P., & van den Bosch, F. G. 2005, 

Mon. Not. R. Astron. Soc, 358, 217 
York, D. G., et al. 2000, Astron. J., 120, 1579 
Zehavi, L, et al. 2004, preprint I astro-p h/0408569 I 
Zheng, Z., et al. 2004, preprint |astro-ph/ 0408564| 



© 0000 RAS, MNRAS 000, 000-000 



